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Abstract. - Two-particle (2-p) excitations such as spin and charge excitations play a key role 
in high-Tc cuprate superconductors (HTSC). On the basis of a parameter- free theory, which 
extends the Variational Cluster Approach (a recently developed embedded cluster method) to 2-p 
excitations, the magnetic excitations of HTSC are shown to be reproduced for a Hubbard model 
within the relevant strong-coupling regime. In particular, the resonance mode in the underdoped 
regime, its intensity and "hour-glass" dispersion are in good overall agreement with experiments. 
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Introduction. — Two-particle (2-p) excitations and 
their corresponding magnetic, charge, optical and pair- 
itig susceptibilities are fundamental for obtaining a micro- 
scopic understanding of the high-Tc cuprate superconduc- 
tor (HTSC) physics, complementing single-particle (such 
as angle-resolved photoemission (ARPES), etc.) experi- 
ments. A key example is provided by the magnetic exci- 
tations: When entering the superconducting (SC) state 
in the high-T^. cuprates, the magnetic excitation spec- 
trum is characteristically and markedly modified: a res- 
onant mode emerges with its peak intensity being highest 
around the wave vector q^AF = (tt, tt) characteristic of an- 
tiferromagnetism (AF) in the undoped parent compound 
Its frequency a;ios(qA_F) follows the doping depen- 
dence of Tc. Away from c^af^ the mode has both a down- 
ward and upward "hour-glass" -like dispersion. A variety 
of experiments in the HTSC, such as photoemission, op- 
tical and tunneling spectroscopies, have been interpreted 
as evidence of interactions of electrons with this mode [7] . 
However its microscopic origin, in particular its role in 
pairing and the more detailed effects arising from the in- 
teractions of charge carriers with this magnetic mode are 
still unclear and intensively debated A prerequisite 

to resolve this debate obviously requires a consistent the- 
oretical description of the neutron resonance mode and. 



more generally, the magnetic excitation spectrum |18j and 
at the same time of the phase diagram, containing the 
competing AF and SC phases. 

In this letter, on the basis of a microscopic theory for 2- 
p excitations, we provide such a consistent description for 
the experimentally relevant regime of the two-dimensional 
(2D) Hubbard model. The essential new points here are 
that our theory for 2-p excitations (e.g. the dynamic spin- 
susceptibility) is (i) parameter-free (given fixed, widely- 
accepted values for the Hubbard model parameters) and 
(ii) is working in the relevant strong correlation regime 
of the underlying Hubbard model. Previous descriptions 
of the magnetic resonance have been obtained by weak- 
coupling |21| and/or semiphenomenological approaches 
|23[|28| reproducing the experimental behavior with ad- 
justable parameters. An accurate description of the infi- 
nite lattice is crucial in order to obtain the magnetic 
resonance which may be considered as a "fingerprint" of 
the AF order in the SC state. Only then are we able to 
differentiate between the competing AF and SC orders in 
the phase diagram. Therefore, the infinite-lattice limit 
(not to be confused with the limit of infinite dimensions) 
has also to be embedded in a controlled description of the 
corresponding susceptibilities. 

We extend the original idea of the Variational Cluster 
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Approach (VGA), which is to extrapolate cluster results 
to the infinite lattice, to the treatment of 2-p excita- 
tions. In our novel approach, the 2-p vertex extracted from 
the corresponding cluster susceptibilities is used to obtain 
the susceptibilities for an infinite lattice. The VGA 
was recently applied to calculate the zero-tempcraturc 
(T=0) phase diagram as well as singlc-particlc excitations 
[30, 32, 33j of the single-band Hubbard model. These re- 
sults succesfuUy reproduced salient experimental features 
such as the electron-hole asymmetry in the doping depen- 
dence of AF and SG phases [31[53] in the HTSG materi- 
als. Also the VGA single-particle excitations were found 
to reproduce characteristic features observed in ARPES 
experiments. In particular, in Ref. [37], the magnitude 
and doping dependence of the SG gap near the nodal and 
antinodal regions was studied in detail, and shown to re- 
produce qualitatively the much-discussed presence of a gap 
dichotomy of the nodal and antinodal SG gaps. Gom- 
bined with the new results for 2-p magnetic excitations, 
presented in this work, a consistent picture emerges, which 
lends substantial support to Hubbard-model descriptions 
of high-Tj. cuprate superconductivity. 

For the appropriate strongly correlated regime {U ~ 8t) 
of the underdoped Hubbard model the resonance is ob- 
tained in a parameter-free calculation and verified to be a 
spin S ~ 1 excitonic bound state, which appears in the SG- 
induced gap in the spectrum of electron-hole spin-flip (i.e. 
S = I) excitations. This will be detailed in our results, 
where we find the doping dependence of (^tcs{(Iaf), the 
energy-integrated spectral weight evaluated at qAF and 
the difference of the magnetic susceptibilities in the SG 
and the normal (N) states to be in qualitative accord with 
neutron scattering data for underdoped YBa2Gu306+x 
(YBGO), where the mode was studied in great detail [l]|6]. 
A spin excitonic bound state has previously been sug- 
gested on the basis of an itinerant picture, most frequently 
invoking a weakly correlated RPA-like form of the dynamic 
spin susceptibility (for a recent reference see, for example, 
Ref. Hi]). As a weak-coupling form it leads to a Fermi- 
liquid like x((7, cj), which is in contrast to some of the 
anomalous dynamics found in neutron scattering experi- 
ments [Hd]. On the other hand, when the 2-p interaction 
and the SG gap arc used as adjustable parameters, it qual- 
itatively accounts for the mode behavior near optimal and 
overdoped regimes |21j . 

Model and Method. — We start from the two- 
dimensional Hubbard model: 

H = -Y^ Ujcl^Cj^ + UJ2 "'t"4 ' (1) 

ija i 

where Uj denote nearest (t) neighbor and next-nearest 
{t' = —0.3t) neighbor hopping matrix elements, c}^ and 
Cjo- are the usual creation and destruction operators, riio- 
their density and U = 8t the local Hubbard repulsion. 
We consider the following matrix expression for the 



transverse spin susceptibility: 

x(Q,*^?„) = x''(Q,*0+ (2) 

with the bosonic Matsubara frequencies — 2rrmT and 
T the temperature. 

Within our embedded cluster approach, the susceptibil- 
ity X depends on the wave vector Q in the reduced Bril- 
louin zone (BZ) associated with the superlattice produced 
by the clusters. In addition, x is considered as a matrix in 
the cluster-site indices i and j and the products in Eq. ^ 
are matrix products. Its definition is given by 

X,,(Q,*0^ rdre-"-^(^--(Q,r)5f(-Q,0)), (3) 
with the mixed representation for the spin operator 

^^"^ R 

Here, the sum is carried out over all clusters, whose posi- 
tion is given by the vector R, Nd is the number of clusters, 
^X+Ti component of the spin operator at site i + R 

and T is the usual imaginary time. 

In Eq. ([2]), an effective particle-hole interaction (ma- 
trix) r is introduced, obtained from an average of the 2-p 
vertex over the additional internal frequencies and mo- 
menta Eqs.([llD as well as the following discourse can 
straightforwardly be generalized to other 2-p susceptibil- 
ities, such as the charge response function, by replacing 
the spin operator in Eq. ([3]) by the corresponding charge 
density operator. 

The "bubble" susceptibility in ^ is obtained as 
a convolution of the "dressed" VGA one-particle (1-p) 
Green's functions, i.e. 

x°-(Q,*a.,''J = --^^ 

(gI^^{K + Q, iu:i + ^^l)Gl<l^^}^, iu:i)+ (5) 
+ F^^^iQ K, ^u;l - tu;i)F*J,'^^{K, *..,{)) . 

^ Notice that, besides this approximation, Eq. ((2]| is the well- 
known Bethe-Salpeter equation, which is in principle exact and valid 
for arbitrary interaction and is not limited to weak coupling as the 
RPA approximation. Neglecting the explicit dependence of the in- 
ternal momenta and frequencies has been found to be a reasonable 
approximation in finitc-T single-cluster QMC calculations 1381 and 
in an extension of the Dynamical Cluster Approximation (DCA) to 
2-p susceptibility calculations 1391 . However, we would like to stress 
that, while the Eq. JJJ appears formally similar to the equation used 
in Ref. 1391 for DCA 2-p quantities, subtle differences occur in build- 
ing up the self-consistent (spin-) response within a cluster: This is 
due to the fact that in the DCA (in contrast to our present VCA- 
type scheme) each cluster site "sees" the same mean-field. We have 
found, as shown e.g. in the "consistency number" a (Eq. Il6l7l l). the 
present scheme to give significantly improved results for comparable 
cluster sizes. 
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Here, iw^ = (2n + 1)ttT denote fermionic Matsubara fre- 
quencies, G^'^^ are normal and F^'^'^ anomalous (i. e. 
characteristic of the SC state) fully interacting Green's 
functions obtained within the VGA and expressed in a 
mixed representation analogous to @, and K is again a 
vector of the reduced BZ. 

Within the VGA [3D], the self-energy of a so-called "ref- 
erence system" is used as an approximation to the one 
of the physical system. The reference system typically 
consists of a subcluster of the original system in which 
single-particle terms can be "optimized" to obtain a sta- 
tionary point of the grand potential for the infinite lattice. 
In a similar spirit, the effective interaction T in Eq. ^ is 
obtained from the corresponding cluster quantity, 

Tiiujl) ^ a [(X°i„,ter(*0)"' - (Xclustcr(*0)~'] (6) 

(cf. also Ref. [10]). Note, that this vertex is defined on the 
cluster and thus independent of Q. In Eq. ([6|), Xciustor 
calculated via the same convolution as in ([5]), but using the 
exact single-particle cluster Green's functions, and Xciustcr 
is the exact cluster susceptibility. The latter quantities are 
obtained via Lanczos exact diagonalization. 

In Eq. ^ , we have introduced a multiplicative constant 
a, which is not a free parameter, but rather it is fixed by 
enforcing the sum rule for the transverse spin susceptibil- 
ity: 

^ XuiQ.^Lol) = {SrS+) (7) 

In our parameter-free approach, the value of a can be used 
to assess the quality of our 2-p scheme: when a is close to 
1, the cluster vertex is a good approximation to the vertex 
of the infinite lattice (see Fig. [TJj). The quality of this 
approximation obviously depends on the cluster size Lc, 
and is expected to become exact for Lc ^ oo, where also 
a — >■ 1. It also certainly depends on the ratio jj and on 
doping. We expect the method to give better results at 
strong coupling, where short-range effects dominate and 
are, thus, well accounted for by modest cluster sizes. Be- 
low, we provide some results in support of this expecta- 
tion. After Fourier transforming over the cluster sites, the 
spin susceptibility x(Q-|-k, Q-f k', ito'^^) acquires a depen- 
dence on two momenta due to the translation symmetry 
breaking introduced by the cluster tiling. The physical 
translation- invariant susceptibility with q = 

Q + k is taken to be its diagonal part x(Q + k,Q + k, iuj^^ ) . 

Results. — Fig. [T^ displays our results for the 
q = Qaf part of the transverse spin spectral function 
Imx^iq = qAF,^) at half-filling, i.e. in the AF phase, 
for different cluster sizes Lc = 2x2, 4x2 and VTO x -s/TO. 
Here, a "finite-size" gap appears , which continuously 
diminishes with increasing cluster sizes (Fig. [T^). It is 
comforting to note that the controlling constant a is very 
close to and also continuously decreasing towards a = 1 
as a function of increasing cluster size (Fig. [TJd). 
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Fig. 1: VGA calculation for the 2D Hubbard model {U = 8t): 
(a) The imaginary part of the susceptibility x*(ciaf,i^) at 
half filling and (b) the sum-rule constant q as a function of 
the cluster size at half-filling; (c) Imx*(q,i^) intensity plot at 
4% hole doping. 

Fig. \Vp shows the corresponding intensity plot for 
lmx^{q,uj) with remnants of the spin-density wave dis- 
persion for a Lc = 2 X 4 "reference" cluster at 4% 
hole doping. Results obtained with different cluster sizes 
reveal that at finite doping "finite-size" effects are of mi- 
nor importance compared to the half-filled situation. We 
attribute this to a screening effect, which renders the 2-p 
vertex F significantly more short-ranged, i.e. more local. 
This means that it can accurately be extracted from the 
exact diagonalization of relatively small clusters. 

Our results for the magnetic response properties in the 
SG state are shown in Fig. [5] Fig. ^ displays a density 
plot of the spin spectral function at x=0.18 doping (3x3 
"reference" cluster). The magnetic resonance emerges in 
the SG-induced gap of = 1 electron-hole (e-h) excita- 
tions when entering this SG doping regime. In Fig. [5^, 
wc plot the intensity in the (w, q)-planc along the diago- 
nal of the 2D-BZ. The celebrated "hour-glass" structure 
observed in the experiments (see, in particular, Ref. [6]) 
is very well reproduced in our calculation. The struc- 
ture has its maximum spectral weight confined to a region 
close to qAF and a dramatic intensity reduction around 
~ 0.8(7r,7r). 

Our results confirm the experimental interpretation put 
forward in Ref. [5]. As shown in a density plot in Fig. [^b, 
we find for the same doping as in Fig. ^ {x = 0.18) a 
typical Fermi surface closed around (tTjTt). Fig. [2j: plots 
the corresponding e-h continuum (i.e. S = I, spin-flip 
c-h excitations), obtained in our VGA calculation. Only 
collective modes below the e-h continuum (red line) can 
actually be detected, because modes within the continuum 
are Landau damped. This e-h continuum corresponds to 
Fig. 4c in Ref. [6]. The continuum threshold exhibits also 
in our case a pronounced minimum in the vicinity of the 
wave vector 2feN — 0.8(7r,7r), which corresponds to scat- 
tering between nodes of the d-wave gap function (Fig. [2b 
gives just one quadrant of the BZ). The minimum in our 
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Fig. 2: (a) Intensity plot of lmx^{q,i^) from the microscopic 
2-p tiieory Eqs. ((2][7|) along q = rj{Tr, ir) displaying the "hour- 
glass" shape, (b) Intensity plot of the low-energy spectral 
weight obtained from the corresponding VGA calculation for 
the 1-p spectral weight displaying the Fermi surface with the 
nodal scattering vector 2fcN — 0.8(7r,7r). (c) Spin-flip electron- 
hole continuum (hatched area: extracted from Eq. Q) with a 
minimum at 2feN. All results are obtained in the SC phase at 
X = 0.18. 



calculation is, however, not so steep as in the idealistic sit- 
uation in Fig. 4c of Rcf. [6], due to correlation effects and 
to a broadening of 0.05t used for the exact diagonalization. 

We would like to emphasize that a similar picture has 
been suggested in RPA-like descriptions of the neutron 
resonance (see, for example, Ref. 121] )■ However, in these 
calculations the d-wave gap amplitude as well as the mag- 
nitude of the effective 2-p interaction have been introduced 
as adjustable parameters in order to reproduce the exper- 
imental energy positions of the resonance mode at (tt, tt) 
and the e-h threshold around 0.8(7r, tt). 

In Figs. [5^ to[2t, we show additional comparisons of our 
calculations with salient features of the neutron scattering 
experiments in underdoped YBa2Cu306-i-x [S]: Fig. [3^ 
displays the difference (ImAx^(q^i;^,a;)) between the spin 
spectral functions in the SC and in the N phases. The 
neutron results by Fong et al. [U] (left panel) are com- 
pared with our results (right panel). Since we are using 
a r = method, our "normal-state" solutions have been 
obtained by forcing the SC Weiss field to be zero in the 
variational procedure. Furthermore, the dopings in the 
experimental and in the theoretical curves do not exactly 
coincide, although both are in the underdoped region. In 
this sense, our comparison is only qualitative. Neverthe- 
less, our calculations reproduce the experimental finding, 
that the enhancement of the spectral weight around the 
resonance peak energy is accompanied by a reduction of 
the spectral weight over a limited energy range both above 
and below ij-[csi<iAF) ■ 

Fig. [3)d compares the energy-integrated spin spectral 
weight at qaf obtained in experiment at various dop- 
ings [H] (left panel) with our theoretical results (right) 
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Fig. 3: Comparisons of our theoretical results with experi- 
ment (reprinted with kind permission from Ref. £l|] "Copy- 
right (2000) by the American Physical Society"): (a) Difference 
between lmx{ciAF , i^) in the SC and normal states (theory: 
X — 0.17); (b) w-integrated spectral weight at qaf; (c) cJi-cs as 
a function of doping. 



in the SC region. The overall doping dependence is simi- 
lar. Finally, in Fig. the doping behavior of the energy 
of the magnetic resonance is compared in the underdoped 
regime (x = 1 corresponds in experiment to optimal dop- 
ing). Again a similar trend is observed. 

Summary and conclusions. In summary, our 
novel theory for 2-p excitations is able to provide an ap- 
propriate description of the resonance mode in HTSC in 
good agreement with experiments. In particular, the cal- 
culated doping dependence of a;ros(qAF), the "hour-glass" 
dispersion of the resonance and its rapid decrease around 
a characteristic wave vector 2fcN, which coincides with the 
distance between nodal points on the Fermi surface, are 
qualitatively consistent with the experiment and support 
the S = 1 magnetic exciton scenario. In contrast to pre- 
vious calculations, our results are obtained in the appro- 
priate strong-correlation regime and contain no adjustable 
parameters. When taken together with earlier results on 
the phase diagram and single-particle excitations, they 
constitute a rather strong support for a Hubbard-model 
description of the HTSC materials. 
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